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We present the results of a numerical search for periodic orbits of three equal masses moving 
in a plane under the influence of Newtonian gravity, with zero angular momentum. A topological 
method is used to classify periodic three-body orbits into families, which fall into four classes, with 
all three previously known families belonging to one class. The classes are defined by the orbits 
geometric and algebraic symmetries. In each class we present a few orbits initial conditions, 15 in 
all; 13 of these correspond to distinct orbits. 
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After Bruns showed that there are 18 degrees-of- 
freedom, but only 10 integrals-of-motion in the dynamics 
of three Newtonian bodies, late in the 19th century, Ref. 
[l| , it has been clear that the three-body problem can not 
be solved in the same sense as the two-bodyone. That 
realization led to Poincare's famous dictum [2j "... what 
makes these (periodic) solutions so precious to us, is that 
they are, so to say, the only opening through which we 
can try to penetrate in a place which, up to now, was sup- 
posed to be inaccessible". Consequently (new) periodic 
three-body solutions have been sought ever since, though 
a significant number were found only after 1975. They 
may be classified in three families: 1) the Lagrange-Euler 
one, dating back to the eighteenth century analytical so- 
lutions, supplemented by one recent orbit due to Moore 
3, 2) the Broucke-Henon-Hadjidemetriou family, dating 
to the mid-1970s [J-Q with periodic rediscoveries of cer- 
tain members of this family 0, and 3) the Figure-8 
family discovered by Moore in 1993, Ref. 0|, rediscov- 
ered in 2000, Ref. [1J| , and extended to the rotating case 
in Refs. hj-nj, (see also Ref. [l8[ and the gallery of 
orbits in Ref. [ill ). 

The aforementioned rediscoveries raise the issue of 
proper identification and classification of periodic three- 
body trajectories. Moore Q used braids drawn out by 
the three particles' trajectories in 2+1 dimensional space- 
time, Ref. [20], to label periodic solutions. This method 
does not associate a periodic orbit with a single braid, 
however, but with the "conjugacy class" of a braid group 
element, i.e., with all cyclic permutations of the strand 
crossings constituting a particular braid. While reason- 
ably effective for the identification of individual orbits, 
braids are less efficient at classifying orbits into families. 

Montgomery [2l[ suggested using the topological prop- 
erties of trajectories on the so-called shape-space sphere 
22] to classify families of three-body orbits. That 
method led Chenciner and Montgomery to their redis- 
covery of the figure-8 orbit [ll[ and informed the present 
study. No solutions belonging to new topological classes 



"higher" than the figure-8 one have been found in New- 
tonian gravity since then, however. 

Here we report the results of our ongoing numeri- 
cal search for periodic collisionless planar solutions with 
zero-angular-momentum in a two-parameter subspace of 
(the full four-dimensional space of) scaled zero-angular 
momentum initial conditions. This subspace is defined as 
that of collinear configurations with one body exactly in 
the middle between the other two, with vanishing angular 
momentum and vanishing time derivative of the hyper- 
radius at the initial time. At first we found around 50 dif- 
ferent regions containing candidates for periodic orbits, 
at return proximity of 10 _1 in the phase space, in this 
section of the initial conditions space. Then, we refined 
these initial conditions to the level of return proximity of 
less then < 10~ 6 by using the gradient descent method. 
Here we present 15 solutions, which can be classified into 
13 topologically distinct families. This is because two 
pairs of initial conditions specify only two independent 
solutions, as the respective members of the pairs are re- 
lated by a simple rescaling of space and time. Before 
describing these orbits and their families we must specify 
the topological classification method more closely. 

Montgomery [2lJ noticed the connection between the 
"fundamental group of a two sphere with three punc- 
tures," i.e., the "free group on two letters" (a, b), and 
the conjugacy classes of the "projective coloured or pure 
braid group" of three strands PB^. Graphically, this 
method amounts to classifying closed curves according 
to their topologies on a sphere with three punctures. A 
stereographic projection of this sphere onto a plane, us- 
ing one of the punctures as the "north pole" effectively 
removes that puncture to infinity, and reduces the prob- 
lem to one of classifying closed curves in a plane with 
two punctures. That leads to the aforementioned free 
group on two letters (a, 6), where (for definiteness) a de- 
notes a clockwise full turn around the right-hand-side 
puncture, and b denotes the counterclockwise full turn 
around the other puncture, see Ref. [llj]. For better 
legibility we denote their inverses by capitalized letters 



= A, b~ 



B. Each family of orbits is associated 
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with the conjugacy class of a free group element. For 
example the conjugacy class of the free group element 
aB contains A(aB)a — Ba. To appreciate the utility of 



FIG. 1: The (translucent) shape-space sphere, with its back side also visible here. Three two-body collision points (bold red 
circles) - punctures in the sphere - lie on the equator, (a) The solid black line encircling the shape sphere twice is the figure-8 
orbit, (b) Class LA butterfly I orbit (I.A.I). Note the two reflection symmetry axes, (c) Class LB moth I orbit (I.B.I) on the 
shape-space sphere. Note the two reflection symmetry axes, (d) Class II. B yarn orbit (II.B.l) on the shape-space sphere. Note 
the single-point reflection symmetry, (e) Class II. C yin-yang I orbit (II. C. 2) on the shape-space sphere. Note the single-point 
reflection symmetry, (f) An illustration of a real space orbit, the "yin-yang II" orbit (II.C.3a). 



this classification one must first identify the two-sphere 
with three punctures with the shape-space sphere and 
the three two-body collision points with the punctures. 

With two three-body Jacobi relative coordinate vec- 
tors, p = ^|(xi -x 2 ), A = ^(xi +x 2 -2x 3 ), there are 
three independent scalar three-body variables, i.e., p ■ A, 
p 2 , and A 2 . Thus the "internal configuration space" of 
the planar three-body problem is three-dimensional. The 
hyper-radius R = y/ p 2 + A 2 defines the overall size of the 
system and removes one of the three linear combinations 
of scalar variables. Thus, one may relate the three scalars 
to a (unit) hyperspace three-vector n with the Cartesian 

components n x — ^r-, n y — 2/ ^* A and n z — The 
domain of these three-body variables is a sphere with unit 
radius J22j, see Ref. [l8[ and Fig. HJa). The equatorial 
circle corresponds to collinear configurations (degener- 
ate triangles) and the three points on it correspond to 
the two-body collisions (these are Montgomery's "punc- 
tures" ) . 

If one disallows collisions in a periodic orbit, then the 



orbit's trajectory on the sphere cannot be continuously 
stretched over any one of these three punctures, and the 
orbit's characteristic conjugacy class is thereby fixed; in 
this sense the topology characterizes the orbit. Thus, pe- 
riodic solutions belonging to a single collisionless family 
are topologically equivalent closed curves on the shape- 
space sphere with three punctures in it. For example, the 
three previously known families of orbits in shape space 
are shown in Ref. [l8| . 

One may divide the orbits into two types according to 
their symmetries in the shape space: (I) those with reflec- 
tion symmetries about two orthogonal axes - the equa- 
tor and the zeroth meridian passing through the "far" 
collision point; and (II) those with a central reflection 
symmetry about a single point - the intersection of the 
equator and the aforementioned zeroth meridian. Sim- 
ilarly, one may divide the orbits according to algebraic 
exchange symmetries of (conjugacy classes of) their free 
group elements: (A) with free group elements that are 
symmetric under a A and b o B, (B) with free group 
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TABLE I: Initial conditions and periods of three-body orbits. xi(0), yi(0) are the first particle's initial velocities in the x and 
y directions, respectively, T is the period. The other two particles' initial conditions are specified by these two parameters, 
as follows, xi(0) = -2:2(0) = -1, x 3 (0) = 0, yi (0) = y 2 (0) = y 3 (0) = 0, i 2 (0) = ii(0), ± 3 (0) = -2i 1 (0), y 2 (0) = yi(0), 
1/3(0) = — 2yi(0). The Newton's gravity coupling constant G is taken as G — 1 and equal masses as mi,2, 3 = 1. All solutions 
have "inversion partners" (mirror images) in all four quadrants, i.e. if ii(0),yi(0) is a solution, so are ±±i(0), ±7/1 (0). Some 
of these partners are exactly identical to the originals, others are identical up to time reversal, and yet others are related to 
the originals by a reflection; we consider all of them to be physically equivalent to the originals. Note that two pairs of initial 
conditions in the same quadrant (II.C.2a and II.C.2b; and II.C.3a and II.C.3b) specify only two independent solutions; see the 
text for explanation. 



Class, number and name x\ (0) yi(0) 



Free group element 



I.A.I butterfly I 0.30689 0.12551 6.2356 

I.A.2 butterfly II 0.39295 0.09758 7.0039 

I.A.3 bumblebee 0.18428 0.58719 63.5345 



I.B.I moth I 
LB. 2 moth II 
LB. 3 butterfly III 
I.B.4 moth III 
LB. 5 goggles 
I.B.6 butterfly IV 
LB. 7 dragonfly 



0.46444 
0.43917 
0.40592 
0.38344 
0.08330 
0.350112 
0.08058 



0.39606 
0.45297 
0.23016 
0.37736 
0.12789 
0.07934 
0.58884 



14.8939 
28.6703 
13.8658 
25.8406 
10.4668 
79.4759 
21.2710 



II.B.l yarn 0.55906 0.34919 55.5018 

II.C.2a yin-yang I 0.51394 0.30474 17.3284 

II.C.2b yin-yang I 0.28270 0.32721 10.9626 

LLC. 3a yin-yang II 0.41682 0.33033 55.7898 

II.C.3b yin-yang II 0.41734 0.31310 54.2076 



abf{ABY 
abf(AB) 2 

b 2 (ABab) 2 A 2 (baBAfba) (B 2 (abABfa 2 (BAbafBA) 



<a(B AB)ab(AB A) 

abAB) 2 A(baBA) 2 B 

ab) 2 (ABA) (ba) 2 (BAB) 

babABA) 2 a(abaBAB) 2 b 

ab) 2 AB B A(ba) 2 B AAB 

(ab) 2 (AB) 2 ) 6 A((ba) 2 (BA) 2 ) 6 B 

b 2 (ABabAB))(a 2 (BAbaBA)) 



babABabaBAf 
ab) 2 (AB A)ba{B AB) 
ab)(ABA)ba(BAB) 
abaBAB) 3 (abaBAbab) (AB Ababf \AB) 2 
abaBABf(abaBAbab){ABAbabf(AB) 2 



elements symmetric under a -f-> b and A -o- B, and (C) 
with free group elements that are not symmetric under 
either of the two symmetries (A) or (B). We have ob- 
served empirically that, for all presently known orbits, 
the algebraic symmetry class (A) always corresponds to 
the geometric class (I), and that the algebraic class (C) 
always corresponds to the geometric class (II), whereas 
the algebraic class (B) may fall into either of the two ge- 
ometric classes. The first examples, to our knowledge, 
of higher topology trajectories on the shape-space sphere 
are the two (new) zero-angular-momentum periodic so- 
lutions reported in Ref. [23j, albeit in a different (the 
so-called Y-string) potential. Here we show only the new 
orbits in Newtonian gravity. 

(I. A) As new members of this class, we present three 
orbits in Table Q] butterflies I & II and the bumblebee. 
We show the butterfly I in Fig. QJb). The butterfly's free 
group element is (ab) 2 (AB) 2 . Note its close relation to 
the figure-8 orbit's free group element (ab)(AB) - both 
orbits belong to this class. We have found two distinct 
butterfly orbits with the same topology (see Table HJ but 
with different periods and sizes of trajectories, both on 
the shape sphere and in real space, see Ref. [18|. This 
kind of multiplicity of solutions is not the first one of 
its kind: there are two (very similar in appearance, yet 
distinct) kinds of figure-8, [14] , 

(I.B) An example of this class of solutions is the moth 
I orbit, shown in Fig. Hfc). We have found a number 
of other solutions that belong to this class of solutions 
with visibly different geometrical patterns on the shape 



sphere and different free group elements; see Table U and 
Ref. [3. 

(II. B) An example of this class of solutions with alge- 
braic symmetry (B), but with only a central geometric 
symmetry, is the yarn orbit (II.B.l), shown in Fig. Hid). 

(II. C) An example of this class without algebraic sym- 
metries is the simplest zero-angular-momentum yin-yang 
I orbit (II. C. 2), shown in Fig. [TJe). There are two dif- 
ferent sets of initial conditions (see Table H]) that lead to 
the same yin-yang orbit in shape space, due to the fact 
that this trajectory crosses the initial configuration on 
the shape-space sphere twice in one period, albeit with 
different velocities. Therefore the two sets of initial con- 
ditions have different energies, so that their periods are 
different, yet both correspond to the same orbit, modulo 
rescaling of the space and time, see Ref. [7]. We have 
found four sets of initial conditions (see Table |T| corre- 
sponding to two distinct (i.e. having different free group 
elements) solutions that belong to this (yin-yang) general 
class. All yin-yang orbits seem to emerge from a single 
quasi-one-dimensional periodic orbit with collisions [18| , 
very much like the Broucke-Hcnon-Hadjidemetriou ones 
emerge from the Schubart (colliding) orbit [24T |. 

In conclusion, we have shown 13 new, distinct equal 
mass, zero-angular-momcntum, planar collision-less peri- 
odic three-body orbits that can be classified in three new 
(and one old) classes. If the figure-8 orbit and its family 
can be used as a benchmark, then we expect each of the 
new orbits to define a family of periodic solutions with 
nonzero angular momentum. We expect our solutions 
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to be either stable or marginally unstable, as otherwise 
they probably would not have been found by the present 
method. 

No three objects with equal masses and zero angu- 
lar momentum, have been found by observational as- 
tronomers, as yet, so our solutions cannot be directly 
compared with observed data [25j . Most of the three- 
body systems identified in observations thus far belong 
either to the (Euler-)Lagrange class, or to the quasi- 
Keplerian Broucke-Hadjidemetriou-Henon class of solu- 
tions. 

Besides obvious questions, such as the study of sta- 
bility, and the search for the associated nonzero-angular- 
momentum solutions, there are other directions for future 
research, such as the nonequal mass solutions [13], the 
general- relativistic extensions of these orbits Ref. [28j . 



as well as the gravitational wave patterns that they gen- 
erate [HH^. 

Having searched in only one two-dimensional section 
of the full four-dimensional space of initial conditions, we 
expect other types of orbits to appear, [single members 
of which have already been seen e.g. the goggles and the 
dragonfly orbits, in Ref. [HI], in different sections of the 
full space of initial conditions. Last, but not least, new 
numerical solutions in the Newtonian potential may lead 
to new analytical solutions, for example of the kind found 
in Ref. [31| after the numerical discovery of the figure-8 
orbit, albeit in the — 1/r 2 potential. Thus, our work may 
shed further light on the three-body problem. 

This work was supported by the Serbian Ministry of 
Education, Science and Technological Development un- 
der Grant No. 01 171037 and No. Ill 41011. 
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